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PREFACE 


A description is presented of the development and preliminary results of the Wave 
Climate Model of the Mid- Atlantic Continental Shelf and Shoreline (Virginian Sea), a joint 
effort of the Virginia Institute of Marine Science (VIMS) and the Langley Research Center 
(LaRC). A review of the shelf geomorphology is also presented, since the most important 
influence on the wave climate of this shelf is the interaction between the ocean waves and 
the various shelf relief elements. 

This Wave Climate Model resulted from extensive modifications made by Victor 
Goldsmith and Joseph M. Colonell of a wave refraction program developed by Dobson 
(Stanford Univ. Tech. Rep. No. 80). The resulting analytical model was adapted to the 
Control Data 6000 computer systems at the Langley Research Center by Vincent R. 

Roland and James R. Schiess. The map projection constructed at the Clark University 
Cartography Laboratory is described by Norman T. Carpenter in appendix A. The 
84 420 depths on this grid were accumulated through the efforts of Carolyn H. Sutton and 
Gaynor Williams. Joseph M. Colonell of the University of Massachusetts, George Grant 
and Maynard M. Nichols of VIMS, and Theodore A. Talay and Andrew R. Wineman of 
LaRC reviewed the manuscript. 

This research was supported by Sea Grant Project No. 5-72 to VIMS from the Office 
of Sea Grant Programs, National Oceanic and Atmospheric Administration, Department of 
Commerce, under Public Law 89-688, and the Commonwealth of Virginia. 
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INTRODUCTION 


As man has increased his usage of the continental shelves of the world, the need to 
understand those physical processes affecting his planning and operations has grown con- 
comitantly. For example, as population pressures have increased shoreline recreational 
demands, the need for design criteria information for shoreline defenses has grown. As 
the need for offshore power plants or port facilities has grown, so has grown the need for 
that environmental information required for site selection and design. In both of the 
cited examples and others, one of the most urgently needed information elements is that 
of surface wave behavior on the continental shelf and fringing coastline. 

In most engineering activities on the coastline and shelf, the designer needs to know 
the temporal and spatial variations in wave energy so that he may select a site or design 
a structure. Experience has shown that the acquisition of sufficiently verified statistical 
information on waves in the open ocean is an enormously expensive and technically diffi- 
cult undertaking. Even if adequate wave information is available for the deep-ocean 
regions, it is necessary to manipulate this information so that the transformations which 
waves experience in transit over the shelf are accounted for. Recent and future theoret- 
ical advances in our understanding of wave generation, coupled with advances in remote 
sensing technology, foster the expectation that adequate deep-sea wave information will, 
in the near future, be routinely observed. The purpose of this report is to supply a 
detailed description of the development of a methodology to investigate the expected 
behavior of waves as they pass over the continental shelf and the resulting wave energy 
distribution along the coastline. This Virginia Institute of Marine Science (VIMS) - Langley 
Research Center (LaRC) Wave Climate Model supplies the interfacing function whereby 
shelf and shoreline wave behavior may be calculated for any specified ocean wave input. 

The model applies linear wave theory and shelf depth information so that a first- 
order approximation of the varying wave path and other parametric changes is obtained. 
Because of the advent of high-speed digital computers such models have become increas- 
ingly common within the last few years (refs. 1 to 21). Basically all these models move 
shoreward a single wave ray at a time across a grid of depths. Wave behavior is deter- 
mined for a number of rays for each specific wave condition (i.e., wave period, height, 
direction, and tide). In fact, Bascom (ref. 22, p. 11) has correctly observed, . . the 
theoreticians have become so bold as a result of the success of their complicated equa- 
tions that there is danger the study of waves will fall entirely into the hands of men who 
have never seen the sea." An attempt is being made to avoid this failing (discussed by 
Bascom) in the present studies through comparisons of the results of the Wave Climate 



Model with vertical aerial photographs of the complex sea surface (e.g., Saco Bay, Maine), 
wide application of the data to real-life problems, and adaptation of this model to future 
developments in space technology (e.g., linking the Wave Climate Model to an Earth 
satellite storm warning system). 

The Virginian Sea Wave Climate Model differs from previous models in the follow- 
ing important elements: 

(1) The model covers a very large geographic area of the continental shelf and 
shoreline, Cape Henlopen, Del., to Cape Hatteras, N.C., an area of 20 000 square n. mi. 
within a single large grid. 1 The importance of this approach is that the resulting graph- 
ical display allows the investigator to visually integrate patterns of wave behavior which 
would escape detection when smaller areas are used; as a result, regional differences in 
behavior within the grid stand out. More detailed studies can then be made on a finer - 
mesh grid in specific subareas by using the wave information from the large grid as 
input to the smaller grid. 

(2) Distortions due to flat representations of the spherical Earth and problems 
resulting from the fact that waves travel great circle paths were overcome by construct- 
ing a Transverse Mercator map projection tangent to the Earth along the center of the 

grid. 

(3) An improved understanding of wave behavior in the area of crossed wave rays is 
now available from the theoretical studies of Chao (ref. 13) and Pierson (ref. 23). These 
studies have been applied to the interpretation of such wave phenomena as curved caustics 
(which occur over continental shelf ridge and swale bathymetry) and straight caustics 
(which occur directly over the margins of the deeply incised channels off the mouths of 
the Delaware and Chesapeake Bays). 

The depth grid used utilized an input of 84 420 depths with a unit cell of 0.5 n. mi. 
on a side. The specified wave input conditions considered nine initial directions for six 
different wave frequencies, two wave heights, and two tidal conditions (for three approach 
directions). In all, 122 separate wave conditions were used with 19 different wave 
parameters computed as output for the whole shelf and adjacent shoreline. 

These wave data, which encompass waves propagated shoreward from deep water 
using a 0.5 n. mi. depth grid, are being used (not reported herein) as input to four 
smaller grids in order to illustrate the usefulness of these data to specific problems. 


*A11 shoreline and shelf geographic locations mentioned in the text are shown on 
the bathymetric map in figure 1. 
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Figure 1.- Bathymetry of Virginian Sea (ref. 24). 
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Figure 2.- Comparison of vertical aerial photograph of Saco Bay, Maine, and wave 
refraction diagram computed for T = 8.0 sec and AZ = 135°. (From ref. 12.) 
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Figure 9.- Three-dimensional computer projection of depth data. 
(Used in preparation of frontispiece.) 
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These grids, consisting of three "second 
order" and one "third order," are described 
schematically in sketch (a). Output from 
the first-order grid is used as input to the 
second-order grid, and so on, the result 
being an increasing density of depths as the 
waves approach shallower depths. 



depth grid 


2d order 
0.25 n. mi. 
depth grid 



depth grid 


Deep 

water 


Sketch (a) 


The three second-order grids are located in the vicinity of Wachapreague Inlet 
(1852 and 1934 bathymetry) on the Eastern Shore of Virginia and adjacent to the area 
from Cape Henry, Va., to the Virginia-North Carolina state line. A third-order grid 
encompasses 6 n. mi. of shoreline along the resort area of Virginia Beach, Va. This 
method allows sufficient detail of wave behavior along relatively short stretches of 
coastline without overburdening the computer with the voluminous depth input informa- 
tion required to propagate the waves shoreward from deep water. 

Still another application of this Wave Climate Model conducted during the present 
study is to delineate the importance of shifting wave patterns on the growth and develop- 
ment of the largely offset coastal inlets along the Eastern Shore of Virginia. Two identi- 
cally sized second-order depth grids have been made up for the area adjacent to 
Wachapreague Inlet. One of these grids contains depths acquired in the 1852 offshore 
survey of the area. The companion grid contains depths acquired in the most recent sur- 
vey (1934). The question then raised (to be discussed later) is, did the wave climate in 
1852, derived by inputting the 1852 bathymetric data, drive the barrier island-inlet sys- 
tem towards its greatly changed present configuration (ref. 25). 

This report presents a description of the VIMS-LaRC Virginian Sea Wave Climate 
Model, a review of the geomorphology of the Mid- Atlantic Continental Shelf, the wave 
computations resulting from 122 wave input conditions, and a preliminary analysis of 
these voluminous data. After detailed analyses of these data, the information will be 
applied to problems which are either currently pressing or soon will be. The computed 
wave information is being summarized by contouring the shelf and shoreline for impor- 
tant computed parameters such as wave height, wave energy, bottom orbital velocity, 
shoreline power gradient, and other parameters for important and commonly occurring 
wave conditions. Shelf areas of "confused seas" are being highlighted and shoreline 
areas of wave energy noted. Thus, several techniques and methods will probably be 
developed and evaluated for using the voluminous output from the Wave Climate Model in 
order to increase its usefulness to managers of continental shelf and shoreline resources. 
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SYMBOLS 


AZ azimuth, deg 

b,/3 distance between adjacent wave rays 

C wave celerity 

d water depth 

g acceleration due to gravity 

H wave height 

K r coefficient of refraction 

K g coefficient of shoaling 

k wave number, 2 tt/L 

L wave length 

T wave period 

of angle between successive wave-front positions (or ray orientations) and 

respective adjacent bottom contours 

MLW mean low water 

SSMO summary of synoptic meteorological observation 

UTM Universal Transverse Mercator 

Subscript: 
o initial 
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ANALYTICAL MODEL 


Development History 

Early adaptations of computerized techniques to the study of waves refracting over 
nonbreaking regions were made by Mehr (ref. 26) and Griswold (ref. 27). Harrison and 
Wilson (ref. 28) expanded Griswold's work into a computer program which produced a 
wave orthogonal trace based upon a grid of wave speed in a study of the area off Virginia 
Beach, Va. Such a grid of wave speeds is awkward because it requires a new grid to be 
computed in the area of interest for each wave period, as well as for each wave direction. 
Also, computations of wave intensities are difficult from a grid of wave speeds. Wilson 
(ref. 2) used a grid of depths which fitted a linear surface computed from four adjacent 
depths to calculate the path of the wave orthogonal. Again wave heights were not calcu- 
lated. Roberts (ref. 29) and Keulegan and Harrison (ref. 9) applied refraction techniques 
to the study of tsunamis in order to explain the intense concentrated destruction from a 
tsunami at Crescent City, Calif. 

Dobson (ref. 4) fitted a quadratic surface to a grid of 12 adjacent depths to calculate 
wave heights. A refraction program based on Dobson's program was adapted to the 
CDC 3600 computer and associated CalComp plotter at the University of Massachusetts 
(refs. 7 and 30). These studies showed a strong correlation between the calculated ocean 
wave refraction patterns, observed wave patterns, and short-term shoreline erosion meas- 
ured in the field over a 3-year interval (ref. 31). The spatial variation in shoreline 
wave energy patterns was determined from wave orthogonal spacing along the shoreline 
and wave height variations for important wave approach directions. However, because 
the effects of bottom friction were not included in the wave intensity calculations, the cal- 
culated wave heights were found to be reasonable in some cases but too large in others. 
This aspect was improved by adopting routines developed by Coleman and Wright (ref. 11) 
based upon equations for calculating bottom friction developed by Putnam and modified by 
Bretschneider and Reid (refs. 32 and 33). Th*s resulting program, with some additional 
improvements, has been adapted to the Control Data (CDC) 6000 computer systems at the 
Langley Research Center and associated software. 

Theoretical Review 

Wave theory applied to wave refraction and its limitations has been succinctly sum- 
marized by Colonell, Farrell, and Goldsmith (ref. 18), from which this section has been 
largely drawn. 

The process by which waves are slowed, shortened, and steepened as they travel 
into progressively shallower water is called shoaling . Typically, shoaling does not occur 
uniformly along a wave front so that, as the wave celerity decreases in accordance with 
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its shorter wave length, the wave front bends as a result of the variations in celerity 
along the front. This combination of shoaling and wave-front bending is called refraction 
and, for purposes of analysis, is regarded as analogous to its optical counterpart in the 
refraction of light rays. With the assumption of small wave steepness (amplitude/wave 
length), linear water wave theory (ref. 34) provides the following expression for the 
celerity of a progressive sinusoidal wave traveling through water; 


C 2=litanh^ 

2 17 L 


( 1 ) 


Rearrangement of equation (1), with C = L/T, gives the following familiar expression 
for wave length: 


L = 


gT 2 

2-rr 


tanh 


2nd 

L 


( 2 ) 


Because the hyperbolic tangent tends toward unity for large arguments (i.e., "deep" 
water), it is apparent that for sufficiently large depths the wave length (and consequently, 
wave celerity) is a function only of the wave period. It is easy to show that deep water 
can be assumed with a maximum computational error of less than 0.4 percent when the 
depth exceeds only one -half wave length. For all conditions except those along the top 
boundaries of the study area, the one -half wave-length criteria would have been sufficient. 
However, because of the shallow shelf in this boundary region, deep water was assumed 
to be greater than one-fourth wave length; thereby, the initial maximum possible compu- 
tation error was increased to 8 percent. 

The effects of shoaling and refraction can be estimated by linear wave theory. For 
example, the propagation of surface waves into shallow water is analyzed by consideration 
of the wave energy between two vertical planes which are orthogonal to the wave crests 
and which intersect with the surface to produce wave rays. Energy is assumed not to be 
transmitted along the wave crest; thus, it is not transmitted across wave rays. If it is 
also assumed that the wave period is constant and that there is no loss or gain of energy 
from reflection, percolation, or bottom friction, then linear wave theory provides the 
well-known result, 


-- K r K s 
H o 

The coefficients of refraction and shoaling are given by 



( 3 ) 


( 4 ) 
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and 


Ks = 


2 cosh kd 
,2kd + sinh 2kd, 


1/2 


(5) 


In this report, the effects of bottom friction are incorporated in H. Koh and Le Mehaute 
(ref. 35) have calculated the transformation of progressive waves from deep water to 
shore using first-, third-, and fifth-order theory, with results differing by only 5 percent. 

Wave refraction diagrams .- As a surface wave travels into shallower water the 
bottom exerts an ever greater influence on it, because of the effect of depth in the deter- 
mination of wave celerity (eq. (1)). For example, at a depth equal to one -eighth wave 
length, the wave is slowed down to 65.6 percent of its deep-water velocity. When the 
crests of a train of waves are not parallel to the bottom contours (lines of constant 
depth), the forward portions of the wave crest decrease in speed in such a way that the 
crests tend to become alined with the bottom contours. (See sketch (b).) Wave refrac- 
tion diagrams are used to determine the way in which a wave of given period will respond 
to the bottom topography. 



The primary objective in wave refraction analysis is to compute the shoaling and 
refraction coefficients (eqs. (4) and (5)) and thus to determine the variation of wave 
heights (eq. (3)) in shoaling water and to deduce therefrom the specific kinematic or 
dynamic property of the waves that is required for any given application. For example, 
wave particle velocities and accelerations are required for wave force computations, 
whereas variations in wave energy along the shoreline are utilized for sediment transport 
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and coastal process investigations. In most cases, such deductions of wave properties 
are made by utilizing the results of linear wave theory. 

According to Wiegel (ref. 36, p. 155), the procedures for the preparation of wave 
refraction diagrams have existed at least since 1937. A good example of such an early 
effort is reference 37. The diagrams are constructed in either of two ways: 

(1) Wave-front method : This method is essentially a map of the successive posi- 
tions of a wave front at given time intervals as it moves shoreward. To determine wave 
height variations along the fronts, it is necessary to construct wave rays (or orthogonals) 
which are everywhere perpendicular to the wave fronts. 

(2) Wave -ray (orthogonal) method : This method is a technique which allows the 
wave rays to be drawn directly, without benefit of the wave fronts, by determination of 
the trajectories of selected points on a wave front. 

Both methods are based on the premise that the water wave refraction phenomenon is 
analogous to optical refraction to the extent that Snell’s law is applicable; consequently, 
it is assumed that 

sin a* Co 

= — ( 6 ) 

sin a i Ci 

where and are the angles between successive wave-front positions (or ray ori- 
entations) and the respective adjacent bottom contours, and wave celerities (function of 
depth and wave period) for these positions are denoted by Ci and C 2 . Preparation of 
a wave refraction diagram by either method is basically a step-by-step application of 
equations (1), (2), and (6) to a wave of prescribed period and initial deep water direction 
over a specified bottom topography. Waves are terminated as they reach the shoreline, 
that is, when the wave height reaches a value of 0.78 water depth. 

The advent of modern digital computational facilities has made manual construction 
of wave refraction diagrams virtually obsolete. Although the wave-front method is the 
usual technique for manual construction, the wave-ray method is generally regarded as 
more amenable to the automated preparation made possible by the digital computer. 

Wave refraction diagrams for the investigations reported herein were computed by means 
of a program based on the Stanford University wave refraction program (ref. 4), which 
was adapted for use on the CDC 6400 and 6600 computer systems and the associated 
CalComp plotter at the Langley Research Center. This program utilizes the wave-ray 
method and requires as input information a grid of bathymetric data for the area of inter- 
est as well as definitive wave characteristics such as height, period, and initial direction. 

Validity and limitations .- Although there is reasonable confidence in the general 
validity of wave refraction analysis procedures just briefly described, it should be empha- 
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sized that the underlying theory has bounds to its validity which are exceeded in most 
routine wave refraction computations. As noted previously, linear wave theory is pred- 
icated with such assumptions as sinusoidal wave forms, small wave steepness, and con- 
stant water depth. Casual observation is sufficient to verify that ocean waves are not 
sinusoidal and that they become steep as they approach the shore. However, the 'Irreg- 
ular" sea surface may be looked upon as a complex combination of sinusoidal waves. 

Also, the assumption of constant depth for application of linear theory is not strictly 
correct. 

Another aspect, the proper interpretation of crossed wave rays (or fronts) in the 
refraction diagrams (i.e., caustics), such as those in many of the wave refraction dia- 
grams, does not appear to be the problem it was once thought to be. Chao, in a thorough 
series of theoretical (refs. 38 and 39), wave tank (refs. 40 and 41), and continental shelf 
(refs. 13 and 42) refraction studies of this caustic phenomena, has made the following 
conclusion for such wave refraction studies (ref. 42, p. 20): "The rays, after escaping 

from the caustic regions, eventually follow the continued ray path and the wave conditions 
are determined by the 0 factor just as if no caustic had occurred except that there has 
been a phase shift, which is unobservable because of the randomness of waves in nature. 
These conditions eliminate the necessity of the evaluations of the waves near a 
caustic . . . ." Although some wave height changes may occur in the waves that pass 
through a caustic region, theoretical and wave tank studies (ref. 41) suggest that such 
changes seaward of the zone of breaking waves may be minimal and well within the bounds 
set by other limiting factors such as depth information. Farrell (ref. 43) was able to 
perform an impromptu field evaluation of the wave refraction computations that he had 
performed in conjunction with his research on the mechanisms which are dominant in 
determining the geomorphology of Saco Bay, Maine. Figure 2 provides a direct compari- 
son of a wave pattern observed in a vertical aerial photograph in Saco Bay and the com- 
puted refraction diagram. The input for the diagram was based on the actual wave con- 
ditions as closely as they could be determined from the photograph. The qualitative 
correlation between photograph and diagram suggests that the computational procedure is 
reasonably valid for this situation which is characterized by a complex shoreline and 
irregular bathymetry. Field measurements of wave directions and heights need to be 
made downwave from the region of caustics in order to verify these and other qualitative 
observations and the theoretical considerations of Chao. One must also recognize that 
this comparison is for a single wave condition at a single site over a relatively small 
geographic region. Higher order wave interactions may possibly cause significant effects 
in the far-field region. However, very little field data on the effects from caustics and 
other wave interaction phenomena are available in the literature. 

Steady-state conditions are assumed for the present model. Calculations are made 
for two tide heights, mean low water and spring high tide, which are assumed constant 
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Y-axis, n. mi. 



Figure 2.- Comparison of vertical aerial photograph of Saco Bay, Maine, and wave 
refraction diagram computed for T = 8.0 sec and AZ = 135°. (From ref. 12.) 



over the region of the model. Constant tide height is a traditional assumption for wave 
refraction models of small geographical regions where wave transient time across the 
region is small compared to a tidal cycle. For wide shallow continental shelves, waves 
approaching at an angle may travel 100 n. mi. in areas with depths less than ^ L; the 
transient time for these longer waves may be as much as 2 hours. In such cases, a 
time-variable tide height may be required to correctly represent the wave climate. Such 
a complication may be unnecessary from a practical viewpoint, however, if tide changes 
cause only small effects on the wave patterns. For example, in most of the Virginian Sea, 
where the mean tide range is approximately 1.1 meters (3.5 feet) the change in sea-level 
elevation would be about 0.3 meter (1 foot) during the maximum time of wave transit 
across the shelf. This study examines the magnitude of change in wave patterns for the 
maximum tidal range of 1.2 meters (4 feet) in an effort to provide guidance for future 
research. 

The present model follows traditional practice in that no winds are assumed over 
the region. For small regions where short wave transient times prevail, wave perturba- 
tions caused by surface winds are small. For larger regions with longer wave transient 
times, "headwinds" may possibly have a significant effect on the dissipation of wave 
height. However, very little is known about the effects of winds on waves; neither the 
effects of winds nor currents are considered in this study. 

The theoretical foundation to combat the inherent shortcomings of linear wave 
theory for refraction computations is provided by Stoker (ref. 44) and an excellent review 
of higher order wave theories by Dean (refs. 45 and 46) and Bretschneider (ref. 32) can 
provide guidance for selection of a more suitable theoretical base for any specific 
problem. 

Despite these and other valid criticisms of conventional wave refraction computa- 
tions, it can be argued that the uncertainties associated with the input data required for 
these computations (i.e., depth and wave information) are large enough that the incorpo- 
ration of such theoretical refinements is generally not warranted for most practical 
applications. However, there remains a need for field verification of computational 
results, especially for problems involving fairly irregular bathymetry. The consider- 
able expense of such an undertaking is of course the major deterrent to its performance. 
Nevertheless, there are occasional qualitative field verifications that serve to inspire 
confidence in the analytical procedures (refs. 31, 43, and 25). 

Clearly, substantial caution must be exercised in the interpretation of wave refrac- 
tion diagrams that are subject to these limitations; however, the limitations themselves 
should not be regarded as implacable obstacles to achieving a better understanding of 
coastal hydraulics through the utilization of wave refraction diagrams. Although these 
diagrams are not data in the usual sense, they can provide information on wave behavior 
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that is a useful complement to actual field investigation. This information is essentially 
qualitative in value but, within the limits imposed by linear wave theory, some quantita- 
tive information is also gained without the expense of costly field collection of wave data. 

Other problems to which a wave analyst must be alerted include the role of Earth 
curvature in determining ray travel directions (ref. 13), the shortcomings of representing 
actual ocean waves by theoretically uniform wave trains of distinct period (refs. 47 
and 48), and the problem of extending waves landward of the breaker zone (refs. 49 to 51). 

Map Projection 

Theoretical development .- The definitive mathematical work on map projections 
was computed by Tissot (ref. 52). This work was applied to standard’ map projections by 
Adams (ref. 53), Deetz and Adams (ref. 54), and Army Map Service (ref. 55). A qualita- 
tive description of map projections can be found in references 56 to 58, and a more math- 
ematical description in reference 59. Application of map projections to problems of 
ocean wave propagation have been discussed by Hardy (ref. 60) and Adamo, Baer, and 
Hosiner (ref. 61). The Transverse Mercator map projection specifically constructed for 
this study used a computer program developed by the U.S. Central Intelligence Agency 
(ref. 62). 

Virginian Sea map projection .- Ocean waves travel great circle paths across the 
spherical Earth. Since wave refraction diagrams are flat representations of portions of 
the round Earth, a certain amount of distortion is normally represented in such presenta- 
tions. Where the areas involved are small and close to the equator such distortions are 
minimal. However, for a large area (e.g., Virginian Sea) such distortions can be quite 
critical. For example, 1° of longitude is approximately lA n . mi. wider at latitude 35° N 
than at 39° N. 2 

Whereas Chao (ref. 38) chose to use spherical coordinates based upon latitude and 
longitude instead of Cartesian coordinates for computing the path of the wave ray, the 
present study employed a specially constructed Transverse Mercator map projection for 
solving the round Earth problem. The projection, which has a square grid superimposed 
on it, was chosen because of the need for a common grid in assembling the voluminous 
but detailed depth information from many sources. 

The Mercator projection takes the form of a cylinder which contacts the surface of 
the Earth along the equator. Latitude, longitude, and surface features are transferred to 
the cylinder in such a way that the angular relationships around each point are preserved 
without deformation. Such a projection is called conformal . The cylinder is then 
unwrapped to form a flat map. The Universal Transverse Mercator has the same char- 
acteristics as the standard Mercator except that a meridian has been made the standard 
line instead of the equator. The extreme areal exaggerations that occur with increasing 
distance from the standard line are avoided by limiting the extent of the standard Mercator 
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projection to 3° of longitude on either side of the standard line and using a separate projec- 
tion for each of 60 zones which cover the whole circumference of the Earth. As a result 
the maximum areal exaggeration is one in a thousand over the whole width of one zone. 

The projection prepared herein differs from the standard UTM in three ways: 

First, to simplify computation it is based on the sphere of equal area rather than the 
spheroid (which is an ellipsoid of rotation and a closer mathematical approximation of 
the somewhat irregular shape of the Earth). Although the diameter of the spheroid is 
smaller at the poles and larger at the equator than that of the sphere, they approximate 
each other quite closely, within 1/4 of 1 percent, in the middle latitude for which the 
Virginian Sea projection was constructed. Second, the nominal scale of this projection 
is true at the central meridian of 75°. Third, the rectangular coordinate system super- 
imposed on the UTM has its origin on the equator 500 000 meters (1 640 420 feet) to the 
west of the central meridian, whereas this coordinate system is in feet with an origin in 
the northwest corner of the area plotted. , Not one of these differences has any effect on 
the use for which this particular projection is intended. 

A detailed description, including the strengths and limitations, of this map projec- 
tion constructed at the Clark University Cartography Laboratory by Norman T. Carpenter 
is given in appendix A. 


Data Input 

Depths.- The most important input to the Wave Climate Model is the depth informa- 
tion. The depth information used in this study was taken directly from original Hydro- 
graphic Sounding Sheets (i.e., boat sheets). These are charts with depths written along 
the survey lines (refs. 63 and 64). Though these depths have all been tidally corrected to 
depth at the time of mean low water, no interpretation or interpolation has been used; that 
is, these are the original surveyed depths. A location map of 60 U.S. National Ocean 
Survey (formerly Coast and Geodetic Survey) Hydrographic Sounding Sheets and other data 
used in this study are given in figure 3. 

Despite the wide usage of these original sounding sheets, few sources of written 
information exist on the accuracy criteria desired and met in these surveys as well as 
the corrections employed or not employed and their justifications. In order to fill this 
critical information gap a study on the accuracy of the depth and navigational positioning 
has been made by Sallenger, Goldsmith, and Sutton (ref. 65). Figure 4 which is from this 
study illustrates the different criteria set by the U.S. Coast Survey and its successor 
agencies for surveys of different dates. The depths at which waves of the periods used 
in this study are first significantly refracted by the sea floor irregularities are plotted 
over the Coast Survey accuracy criteria to give an indication of the depth errors influ- 
encing the wave climate model. Only four of the charts which were used in the depth 
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Idenlifi cation 
in fig. 3(a) 

Nu ruber 

Date 

Scale 

Comments 

1 

Map 

237 

1850 

1:400 000 

Soundings only along lines indicated; soundings in fathoms (MLW) 

2 

674 

1859 

1 : 200 000 

Soundings only along lines indicated; soundings in fathoms (MLW) 

3 

965 

1868 

1:40000 

Soundings appear to be in fathoms; age and numerical system used confirm this 

4 

1721 

1886 

1 : 200 000 

Soundings in fathoms, only along lines Indicated (MLW) 

5 

3774 

1915 

1 : 20 000 

Soundings in feet (MLW) 

6 

4089 

1919 

1 : 40 000 

Soundings in feet (MLW); 76° west, soundings are barely readable 

7 

4193 

1921 

1 : 40 OOO 

Soundings in feet (MLW); soundings decrease by a factor of 2 eastward of 75°30' 

3 

4194 

1921 

1 : 40 OOO 

Soundings in Teel (MLW); discontinued at 75°30' east because area is covered tay more recent map 

9 

4255 

1922 

1 : 120 000 

Depth in fathoms (MLW); soundings diminish toward cast 

10 

4286 

1922 

1 : 40 OOO 

Soundings in feet (MLW) 

11 

4944 

1922 

1 : 120000 

Soundings it) fathoms (MLW) 

12 

4951 

1929 

1 : 20 OOO 

Soundings in feet (MLW); shoreline soundings 

13 

5348 

1933 

1 : 40 OOO 

Soundings in feet (MLW) 

14 

5349 

1933 

1:20000 

Soundings in feet (MLW) 

15 

5350 

1933 

1: 120000 

Depth in fathoms (MLW); soundings become sparse pasl 74°50' 

16 

5351 

1933 

1 : 40 OOO 

Soundings in feel (MLW) 

17 

5352 

1933 

1 : 10 000 

Soundings in feel (MLW| 

13 

5353 

1933 

1: 40 000 

Soundings in feet (MLW) 

19 

5354 

1933 

1 ; 20000 

Soundings in feet (MLW); becomes sparse toward edges 

20 

5355 

1934 

1:40000 

Soundings in feet (MLW) 

21 

5357 

1933 

1:20 000 

Soundings in feet (MLW) 

22 

5673 

1034 

1:40 000 

Soundings in feet (MLW) 

23 

5674 

1934 

1 : 40000 

Soundings in feet (MLW) 

24 

5702 

1934 

1: 40000 

Soundings in feet (MLW) 

20 

5703 

1934 

1:20000 

Soundings in feet (MLW) 

2G 

5704 

1934 

1 : 20000 

Soundings in feet (MLW) 

27 

5713 

1927 

1 : 120 000 

Soundings in fathoms (MLW) 

28 

0710 

1934 

1: 40 000 

Soundings in feet (MLW) 

20 

5770 

1935 

1:40000 

Soundings in feet (MLW) 

30 

5771 

1934 

1: 40000 

Soundings in feet (MLW) 

31 

5988 

1935 

1 : 40 OOO 

Soundings in feet (MLW) 

32 

5989 

1835 

1 : 40 000 

Soundings in feet (MLW) 

33 

5990 

1935 

1 : 80000 

Soundings in feet (MLW) 

34 

5991 

1935 

1 : 40 OOO 

Soundings in feet (MLW) 

35 

5992 

1935 

1 : 40 OOO 

Soundings in feet (MLW) 

36 

5993 

1935 

1 : 40 OOO 

Soundings In feet (MLW) 

37 

0995 

1935 

1 : 120 000 

Soundings in fathoms (MI.W) 

38 

6219 

1937 

1 : 120000 

Soundings in feet (M7.W| 

39 

6344 

1938 

1 : 40 000 

Soundings in feet (MLW) 

40 

6345 

1938 

1 : 80 000 

Soundings in fathoms (MLW) 

41 

6595 

1940 

1 : 40 OOO 

Soundings in feet (MLW) 

41 

6976 

1947 

1:40 000 

Soundings in feet (MLW) 

42 

7703 

1948 

1: 10000 

Soundings in feet (MLW) 

43 

8218 

1954 

1:25 000 

Soundings In feet (MLW) 

44 

8711 

1962 

1: 10000 

Soundings In feet (MLW) 

45 

8764 

1927 

1:40000 

Soundings in feet (MLW) 

46 

8809 

1963 

1:20 000 

Soundings In feet (MLW) 

47 

8810 

1963 

1 : 40 000 

Soundings in fathoms (MLW) 

48 

9137 

1970 

1:40 000 

Soundings In feet at 2-foot intervals 

49 

9155 

1970 

1:40 000 

Soundings in feet at 2 -foot intervals 

50 

9171 

1970 

1:40000 

Soundings in feet at 2-foot intervals 

51 

9231(0) 

1971 

1 : 80 000 

Soundings in Fathoms (MLW) 

52 

9231(C) 

1971 

1 : 80 000 

Soundings in fathoms (MLW) 

53 

9243 

197 J 

1 : 80 000 

Soundings in fathoms (MLW) 

54 

Bathymetric Map 
O807N-55 

1967 

1 : 125 OOO 

Soundings in fathoms (MLW) 

55 

0807N-56 

1967 

1: 125 000 

Soundings in fathoms (MLW) 

56 

0807N-57 

1967 

1: 125000 

Soundings in fathoms (MLW) 

57 

Archive 

715025.007 

1971 

1 : 40 000 

Soundings in feet (MLW) 

58 

715025.008 

1971 

1 : 40 000 

Soundings in feet (MLW) 

59 

715025.009 

1971 

1 : 40 000 

Soundings in feet (MLW) 

60 

71504.010 

1971 

1 : 40 000 

Soundings in feel (MI.W) 


(b) Index. 

Figure 3.- Concluded. 
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accumulation were surveyed prior to 1915, and only three of these charts were surveyed 
prior to 1870. These charts (prior to 1915) used for the model were surveyed where the 
depths did not exceed 85 feet. 


tt rn T, sec 



Depth 

(a) Error criteria chart. 


Year 

Criteria 

Circa 1844 

None 

Circa I860 

"The allowable error at .sounding-line crossings was not to be more than 3 percent of the depth, with a 
limiting error of 5 percent." (See ref. 64, p. 218.) 

1878/18B3 

"Lines of soundings at their crossings were not to exceed 'in depths of 15 leet and under, two-tenths of a 
toot; between depths of 15 and 30 feet, three-tenths; 30 and 48 feet, five-tenths; between 48 and 72 feet, 
three-fourths of a toot; between 72 and 96, one foot and a half; and between 96 and 150 feet, two feet. In 
the sea-depths the limit of error should not exceed 1 percent . . . (from instructions to surveyors 

1878 and continued in 1883, as given in ref. 64, p. 221). 

1894 

", . . the admissible percent of error at sounding-line crossings was a maximum of 1.5 percent of the 
depth at that point." (See ref. 64, p. 224.) 

1957 

"Maximum errors: (1) 0 to 11 fm. (0 to 20 m.): 1.0 ft. (0.3 m.); (2) 11 to 55 fm. (20 to 100 m.): 0.5 fm. 
(1.0 m.); (3) 55 fm. (100 m.) and deeper: one percent of depth.” (See ref. 63, p. 20.) 


(b) Accuracy criteria. 

Figure A.- Sounding error criteria for NOS (C&GS) Sounding Sheets from references 64 
and 63 as compiled in reference 65. Depths at 1/4 wave length for wave periods 
used in this study are also shown. 
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Approximately 100 000 of these original uninterpolated depths were transferred 
from the 61 sounding sheets and other data, using latitude and longitude, onto a Trans- 
verse Mercator map projection 2.4 by 1.2 meters (8 by 4 feet), specially constructed for 
the present study. 2 Then 84 420 of these depths were read from the map grid in lines at 
0.5 n. mi. intervals spaced 0.5 n. mi. apart and punched on cards. 

In order to provide a check on this depth accumulation process, a computer pro- 
gram was written which plotted 420 east-west bathymetric profiles. Seven typical pro- 
files are illustrated in figure 5. By adjusting the spacing between the profiles it was 



Figure 5.- Seven of the 420 computer-plotted east-west bathymetric profiles 
selected at intervals of 30 minutes of latitude. 


possible to compare both adjacent data points and adjacent profiles. Thus, mispunched 
numbers, computer cards out of place, and incorrectly read depths, which may be 
expected to occur within a large 84 420 depth array, could be visually located. Correc- 
tions were then made on a point -by -point basis by reference to the original Hydrographic 
Sounding Sheets. This comparison was repeated until all the depths on the tape were 
determined to be correct. 


2 An addendum to this study, a detailed map of the Virginian Sea, was made from 
these accumulated depths which were contoured at 3. 7- meter (12 foot) intervals out to 
150 meters (492 feet) and at 30.5-meter (100 foot) intervals out to 305 meters 
(1000 feet). (See ref. 24.) A black-and-white copy of this 1.22- by 0.61-meter (4 by 
2 foot) multicolored map, with the location names used in this report, is shown as 
figure 1. 
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Once the data had been checked, they were input on punch cards, but the number of 
errors that occurred due to the necessary handling of these 5628 cards indicated that a 
more reliable system was necessary. The data were then transferred to magnetic tape, 
but even the tapes were subject to handling errors. The system finally chosen was to 
record the bathymetry on a data cell; this proved to be the most reliable of the three 
techniques. 

An additional test was performed to determine whether the density of depths (i.e., 
0.5 n. mi. spacing) was sufficient to accurately portray a region of this bathymetric com- 
plexity with respect to the regional refraction patterns. It was determined that, by using 
an existing grid (ref. 31) on a much smaller region, very little difference would exist in 
refraction patterns between this depth density and one of approximately half the spacing 
(which would require approximately four times as many depths). (However, since more 
detailed studies in shallow depths and along the shore require an increase in depth den- 
sity, additional small grids are being used in areas of such studies. The output from the 
larger, regional, first-order grid is used as input at the seaward intersection with the 
smaller grids.) 

The equations used in propagating the wave ray across the grid of depths is 
described by Dobson (ref. 4, pp. 14 to 16) in terms of Cartesian coordinates and time. 
Essentially, the ray is propagated forward as a function of the celerity, which in turn is 
related to the depth coordinates. Start from a known point in a given initial direction and 
propagate the ray forward in a series of small steps by using a complex iteration proce- 
dure that keeps the ray curvature within arbitrary set bounds. If the curvature falls out- 
side these bounds, such as in areas of complex high-relief bathymetry, the ray is termi- 
nated and a message printed. Since the likelihood of many of these small steps falling 
exactly on the grid points is small, a depth interpolation procedure is required. Dobson 
(ref. 4) used a second- degree polynomial employing the method of least squares to 
describe the local surface of fit for 12 adjacent depths. This method is an improvement 
over a linear surface fit (computed by Wilson (ref. 2), May and Tanner (ref. 20), and 
others) in that the ray separation factor b can be calculated; b is needed for deter- 
mining the wave intensity. A statistical test is run on the fit of the quadratic surface to 
the local depth configuration. An unrelated test of the ability of such surfaces to accu- 
rately portray the original configuration determined that 12 depths was the optimum num- 
ber and configuration to use in such a quadratic surface. (See ref. 66.) 

These complex mathematical routines are only as good as the input depth informa- 
tion and there are still numerous U.S. coastal areas (including much of the depth data 
used in this study) whose latest depth information predates electronic navigation and pre- 
cision depth records. (See ref. 65.) 
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Wave conditions .- The second major input to the Wave Climate Model was a wide 
variety of wave conditions. Approximately 200 to 250 wave orthogonals were propagated 
shoreward from deep water 3 for each of 122 wave conditions. The wide variety of wave 
conditions, shown in table I, was chosen in order to model as many different combinations 
as possible of wave period, direction, height, and tidal conditions from amongst the infi- 
nite variety of conditions that occur in real life. 

TABLE L- INITIAL DEEP-WATER WAVE CONDITIONS 


Tide 

Wave direction, 
deg 

Wave periods, 
sec 

Wave height 

m 

ft 

m 

ft 



0 

4,6,8,10 

0.61,1.83 

2,6 



22.5 

4,6,8,10,12,14 

.61,1.83 

2,6 



45 NE 

4,6,8, 10, 12, 14 

.61, 1.83 

2,6 



67.5 

4,6,8,10,12,14 

.61, 1.83 

2,6 



90 E 

4,6,8,10,12,14 

.61,1.83 

2,6 



112.5 

4,6,8, 10, 12,14 

.61, 1.83 

2,6 



135 SE 

4,6,8, 10,12, 14 

.61,1.83 

2,6 



167.5 

4,6,8, 10, 12, 14 

.61, 1.83 

2,6 



180 S 

4,6,8, 10, 12, 14 

.61, 1.83 

2,6 



45 

4,6,8,10,12,14 

.61 

2 



90 

4,6,8, 10, 12, 14 

.61 

2 



135 

4,6,8, 10, 12, 14 

.61 

2 


Thus, a "library file" of a wide variety of wave conditions is accumulated that can 
be used in conjunction with other geological, biological, and chemical studies of the shelf 
and shoreline and as an aid to resource managers charged with choosing sites for off- 
shore ports and shoreline defense structures. 

One might well ask the question as to why this "scatter-gun” approach with respect 
to wave input conditions. Why not zero in on just the most significant waves for calcu- 
lating the wave parameters? There are two reasons why a wide variety of wave condi- 
tions is calculated. First, anyone can testify to the almost infinite variety of wave con- 
ditions that may occur over a long span of time. Second, data for determining the 
precise percentage of time that a given wave condition will occur are presently unavaila- 
ble in most areas. The large frequency of conditions is also needed in order to calculate 

3 Waves propagating landward from deep water will be slowed to 0.996C at 
1 i o .. 

d = - Lo. Similarly, at d = L 0 , the waves are traveling at 0.92C o , and at d = - L 0 , 

they are traveling at 0.66C o . Therefore, the slowing down of the wave is a gradual 
process, and for the purpose of starting the 14- second waves, deep water was considered 

to be where the waves were not appreciably affected by the bottom or at depths = ~ L 0 . 
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parameters such as the mean wave height at a shelf location and the total shoreline wave 
energy along a stretch of coast during an average year. This could be easily calculated 
by summing up, based on frequency of occurrence, the data for a given location from each 
of the calculated wave parameters. Also, in order to determine the effects of storm 
waves along a shoreline, the sequence of weather fronts and resulting storm -generated 
waves is needed (ref. 67). 

Such data come from one of four sources: (1) wave measurements by instruments 
in deep water, (2) wave measurements by instruments in shallow water and along the 
shoreline (i.e., on piers, anchored buoys, etc.), (3) ship wave observations compiled by 
U.S. Naval Oceanographic Office by 10° squares called Marsden squares, and (4) wave 
hindcast calculations. None of these methods has produced data considered adequate for 
this area. 

Data from source (1) are quite rare and, where available, are generally of insuffi- 
cient duration to be statistically valid. Summaries of shallow water wave measurements 
calculated from coastal wave gages were found to bear little relationship to individual 
shipboard wave height and period observations off the east coast of the United States 
(ref. 68, p. n-679), as shown in figure 6. These authors further concluded that if ade- 
quate data were available from shipboard wave observations, wave refraction methods 
would be useful in determining shallow water and shoreline wave parameters. Further- 
more, no procedures presently exist for propagating waves seaward, which use wave 
parameters determined from coastal gages as input. 



Figure 6.- Comparison of cumulative shipboard wave height distribution (Marsden square) 
with Atlantic City, N.J., gage data. Numbers in parentheses are number of observa- 
tions in that subsquare (SS) . (From ref. 68.) 

Ship wave observation data are not the final answer either for determining the per- 
centage frequency of occurrence for a given wave condition, as there are several inherent 
biases built into the present data collection system. Several of these biases, such as the 
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awkward computer forms, are discussed by Harris (ref. 69). Another bias is suspected 
from the interpretation of summary graphs of ship wave observations from the Marsden 
squares adjacent to Cape Cod, Mass. (ref. 31) and the southern portion of the Virginian 
Sea (fig. 7). In these two area summaries, the dominant waves on an annual basis appear 
to approach from the west, despite the proximity of land to the west and more than 
3000 n. mi. of ocean to the east. One possible explanation for this suspected bias is 
related to the fact that the ship wave observations are recorded as part of a voluntary 
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program; ships tend to avoid extreme wave conditions, and when they do encounter severe 
conditions, the assigned observer might find that he has more important duties to perform 
than filling out wave forms. Nevertheless, ship wave observations appear to be the best 
information available at present for summing up individual wave conditions. These data 
have been used with some success in making littoral drift calculations along the coast of 
Florida (ref. 70). 



(b) Direction as a function of height. 
Figure 7.- Concluded. 
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The final method used in summarizing wave conditions is wave hindcast calcula- 
tions. All these methods are discussed elsewhere (refs. 32 and 71). Hindcast calcula- 
tions using the Bretschneider-revised Sverdrup-Munk {significant wave) method has been 
computed for four stations along the U.S. East Coast, including one adjacent to the 
Chesapeake Bay entrance (fig. 8), by using data from weather stations for the 3 -year 
period, 1948-1950 (ref. 72). Note the large discrepancy between the ship wave observa- 
tions and this wave hindcast data (figs. 7(b) and 8). Important considerations in this dis- 
crepancy are the two major assumptions used in wave hindcasting: (1) deep water for 
360° around the hindcast station and (2) the meteorological conditions in the 3-year 
period, 1948-1950, are representative of long-term weather conditions. 



Figure 8.- Deep-water wave hindcast data adjacent to Chesapeake Bay entrance. 
(From ref. 72.) Wave rose showing percent of time waves of different 
heights occur from each direction. 


In order to determine the range of parameters used in this study, the wave condi- 
tions in table I were checked against the hindcast and Marsden square data and were 
found to represent 99.11 percent of the hindcasted wave conditions calculated for the 
Virginian Sea area by Saville in reference 72 (table H) and greater than 59.4 percent of the 
ship wave observations reported in the Marsden square summary for Marsden subsquare 
116-55 which covers the Virginian Sea area (table HI). Most of the remaining observed 
wave conditions (29.2 percent) were from the west; therefore they were not deemed sig- 
nificant for this study because of the limiting fetch and other assumptions in this study. 
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TABLE II.- HINDCAST WAVE DATA FOR OCEAN ADJACENT TO 
CHESAPEAKE BAY ENTRANCE 


^3-year period, 1948-1950; from reference 72] 


AZ 


Percent frequency of occurrence at - 



T = 4, 5 sec 

T = 6, 7 sec 

T = 8, 9 sec 

T = 10, 11 sec 

T = 12, 13 sec 


Total 

N 

NNE 

0.03 

.182 

0.334 

0.182 

0.046 

0.06 


0.03 

.804 

NE 

.958 

1.81 

5.81 

7.0 

3.85 

3.25 

22.678 

ENE 

1.25 

5.61 

11.37 

12.6 

7.2 

5.26 

43.29 

E 

1.49 

5.58 

5.35 

3.52 

1.32 

1.26 

18.52 

ESE 

.87 

2.13 

1.11 

.62 

.56 

.55 

5.84 

SE 

.44 

1.32 

.59 

.53 

.55 

.85 

4.28 

SSE 

.32 

.85 

.30 

.17 

.24 

.21 

2.09 

S 

.45 

.64 

.09 

.14 

.03 

.23 

1.58 

Totals 

5.99 

18.274 

24.802 

24.626 

13.81 

11.61 

99.112 


Direction 


T = <5 sec 


8.2 

5 


TABLE III.- SHIP WAVE OBSERVATIONS 
[Marsden Subsquare 116-55] 


Percent frequency of occurrence at - 


5 to 7 seclT = 7 to 9 seclT = 9 to 11 sec T = 11 to 13 sec 


13 to 15 seel 


1.3 


Total calm . . 
Total west . . 
Total unknown 


Total 


99.8 
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Data Output 

For each of the 122 wave conditions (e.g., combinations of wave directions, periods, 
heights, and lunar and storm tides), two kinds of computer output are generated: 

(1) The wave refraction diagrams show the directions of the waves as they are 
propagated landward under specific wave input and depth conditions. Instead of wave 
fronts, these diagrams display wave rays (i.e., orthogonals) which are everywhere per- 
pendicular to the wave fronts. Seventy of these wave diagrams are presented in appen- 
dix B (figs. B1 to B70). 

(2) The computer printout gives the computations of different wave parameters at 
intervals of 609.6 to 914.4 meters (2000 to 3000 feet) along each wave ray. The rays 
were input in deep water at a spacing of 1 n. mi.; approximately 200 to 250 wave rays 
were used in the wave computations for each condition. The parameters printed for each 
of the conditions in table I are 

(1) Angle of wave approach 

(2) Depth 

(3) Wave length 

(4) Wave celerity 

(5) Refraction coefficient 

(6) Shoaling coefficient 

(7) Wave height (including effects of bottom friction) 

(8) Maximum horizontal component of wave orbital velocity at the bottom 

(9) Maximum horizontal component of wave orbital velocity at middepth 

(10) Total wave energy 

(11) Longshore wave energy 

(12) Wave steepness 

(13) Group velocity 

(14) Total wave power 

(15) Longshore wave power 

and the parameters printed when the wave reaches shore are 

(16) Total wave power gradient 

(17) Longshore wave power gradient 

(18) Longshore drift current velocity 

(19) Longshore sediment drift (five different equations) 
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Parameters (1) to (7) were computed based on reference 4, parameters (8) to (17) 
were from reference 36, and parameter (18) is from reference 34. These computations 
may generate as many as 50 000 values for each wave parameter along each ray for each 
wave condition. These data may be used to contour the sea floor in order to delineate 
areas of high wave-induced current motion and wave energy as well as areas of "confused 
seas." Parameters (8) and (9) are used in computing the drag and inertial components of 
wave forces, and parameters (16) and (17) are used in delineating areas of wave energy 
concentration (i.e., areas which may experience larger rates of wave scour or beach 
erosion). 

Longshore sediment drift was calculated by five different equations taken from ref- 
erences 73 to 78. There was an extremely large scatter in longshore drift values as cal- 
culated by the different equations, with most results appearing to be many times too 
large. These results support the widely recognized need for a significantly greater 
understanding of longshore processes. 

Special Computer Requirements 

The application of the refraction model to such a large geographic area is unique 
and was made possible by the capacity to store a large amount of data. Storage require- 
ments for this study were 330 000 octal locations. Of these, approximately 300 000 loca- 
tions were required for the bathymetric data, and the remaining 30 000 locations were 
required for the program, plotting, and tape generation routines. 

The data were run in sets with a set consisting of one tidal height, one wave direc- 
tion, and six wave periods. The approximate computer requirements for a set were 


Total computer storage field length (octal) 330 000 

Computer time for the CDC 6600 computer, sec 1 250 

Operating system calls 47 000 

LaRC computer system cost $650 


A total of 21 sets were run. A total cost including program adaptation and checking was 
approximately $30 000. 

A number of subsequent modifications have been made to the computer model but 
were not incorporated at the time this study was made. Among these was a method which 
can reduce the computer storage requirements to one-fourth of that presently required. 
This method is called the random-access method, and a peripheral disk is used to store 
the data in modules which can then be retrieved individually; thus, data not required are 
bypassed (ref. 79). This method also reduces cost by approximately 65 percent and 
should allow greater use of the program with computers of more limited storage capacity. 
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SHELF GEOMORPHOLOGY - A REVIEW 


Ever since Uchupi's detailed studies of the continental shelf along the east coast of 
the United States (ref. 80), an increasing number of studies have focused attention on the 
various shelf relief elements. Much of this work is discussed in reference 81. Refer- 
ences 82 to 117 are examples of studies which are of direct interest to the Virginian Sea 
area. Also, references 118 to 139 are studies of significance to the origin of the shelf 
geomorphology of the Virginian Sea area. A classic study on the shoreline effects result- 
ing from offshore changes is Jordaan (ref t 140). 

Many of these studies are aimed primarily at shedding light on the controversy 
concerning the origin of the shelf relief elements - that is, are these features relict or 
presently hydraulically active, or a combination? If relict, how much have they been 
modified? It is not the direct purpose of this study to explore this controversy, but 
merely to point out the importance of complex interactions between the ocean waves pass- 
ing over the shelf and many of these shelf relief elements, which result in a nonuniform 
wave energy distribution over the shelf and along the shoreline. This relationship is dis- 
cussed in detail in the section "Study Results." In this section, the purpose is merely to 
elaborate on those geomorphic features which are significant to the wave climate of the 
Mid- Atlantic Continental Shelf and Shoreline. These features are shown in the bathymet- 
ric map of the Virginian Sea (fig. 1) and also in figure 9, a three-dimensional computer 
projection of the depth data. 

Seven east-west bathymetric profiles at intervals of 30 minutes of latitude taken 
from the 0.5 n. mi. depth grid is shown in figure 5. Two important aspects of the profiles 
for this study are the great width and relatively shallow nature of this portion of the con- 
tinental shelf. The abrupt increase in gradient at the shelf edge is between depths of 61 
and 91 meters (200 and 300 feet) and is located as much as 60 n. mi. from shore. The 
distance from shore at which the ocean waves of different period begin to be appreciably 
affected by the sea floor is shown in figure 10. Thus, a great expanse of the continental 
shelf, and superimposed relief elements, is available for influencing ocean wave behavior. 

A closer examination of these profiles (fig. 5) and the detailed bathymetric map of 
the sea floor (fig. 1) reveals that the shelf surface is not a smooth plain but instead con- 
sists of numerous irregularities. These irregularities may be divided into two groups: 

(1) Large-scale morphogeometry consists mainly of erosional forms cut into the 

shelf such as terraces, channels and valleys, and shelf-edge canyons 

(2) Small-scale shelf relief elements consist of low relief features (i.e., less than 

9.144 meters (30 feet)) of probable depositional origin, most notably ridge 
and swale bathymetry and arcuate (e.g., cape- associated) shoals 
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Figure 9.- Three-dimensional computer projection of depth data. 
(Used in preparation of frontispiece.) 
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Figure 10.- Depths at which waves of different periods begin to be appreciably 

affected by Virginian Sea Floor. 

Whereas the origin of group (1) features is directly related to a lowered sea level, 
group (2) features probably formed since the last rise in sea level under the present shelf 
hydraulic conditions. The most recent eustatic sea level lowering reached its maximum 
extent approximately 15 000 years ago on the Atlantic Continental Shelf. Eustatic sea 
level has been within 1.8 meters (6 feet) of its present level approximately 30 000 to 
35 000 years ago and for the last 4000 years (ref. 128). However, tectonic events may 
have severely altered this sequence of sea level changes in this area. (See refs. 141 
to 143.) 


Large-Scale Morphogeometry 

Terraces .- The depths of the outer edge of the prominent shelf terraces determined 
from an east- west profile along 37° latitude from the mouth of Chesapeake Bay out across 
the shelf to Norfolk Canyon (from ref. 24) are given in table IV and depths of these ter- 
races are compared with the depths of other prominent terraces along the East Coast shelf 
given in reference 80. The most pronounced terraces adjacent to Chesapeake Bay are at 
24, 30, 40, and 86 meters (78, 100, 132, and 282 feet). 

The presence of these terraces on the sea floor indicate a step -like bathymetric 
profile. The effect of the steeper portions of the profiles on the incoming waves will 
depend primarily oh the angle of wave approach to these rises. However, even the steep- 
est rises have relatively low-gradient slopes. The slope is 0°07 r 19'’ for the rise between 
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TABLE IV.- DEPTH TO OUTER EDGE OF TERRACES ON 
THE CONTINENTAL SHELF AND SLOPE 


Depth of outer edge of terraces, m (ft), at - 

Chesapeake Bay, 
Va. 

(a) 

Martha's Vineyard, 
Mass. 

(b) 

Atlantic City. 
N.J. 

(b) 

Onslow Bay, 
N.C. 

(b) 

Savannah, 

Ga. 

(b) 

Cape Kennedy, 
Fla. 

(b) 

Miami , 
Fla. 
(b) 






— 

10 (33) 
15 (49) 
18 (59) 

7 to 18 (24 to 60) 
16 (54) 

24 (78) 

30 (100) 







20 (66) 

20 (66) 

20 (66) 
25 (82) 
30 (98) 


25 (82) 



30 (98) 


35 (115) 


33 (108) 
40 (131) 
45 (148) 
50 (164) 
55 (180) 

33 (108) 


40 (131) 

40 (131) 
45 (148) 

40 (131) 

40 (132) 

43 (141) 



50 (164) 


57 (188) 

55 (180) 
63 (207) 





63 (207) 


62 (203) 



67 (220) 

65 (213) 




70 (230) 
80 (262) 

70 (230) 
80 (262) 

86 (282) 

80 (262) 

83 (272) 
95 (312) 

80 (262) 


106 (348) 


100 (328) 
120 (394) 




120 (394) 





125 (410) 





130 (426) 
140 (459) 











158 (518) 

170 (558) 










175 (574) 




183 (600) 
244 (800) 


200 (656) 




210 (689) 











a Reference 24. 
Reference 80, table 2. 


depths of 87.8 and 62.2 meters (288 and 204 feet) as compared with a slope of 0°01’58 M 
for the total shelf landward of the depth contour of 62.2 meters (204 feet). (See ref. 24.) 

Subaqueous stream drainage . - Generally oriented perpendicular to the strike of the 
terraces, the major relief features remaining from the Pleistocene stream drainage are 
the shelf valleys at the mouths of Delaware and Chesapeake Bays. (However, Swift 
(ref. 115) has suggested that the Delaware shelf valley is an estuary retreat path and not 
a drowned river valley.) Both these southeast-oriented valleys have a pronounced influ- 
ence on the wave refraction patterns, with areas of confused seas forming over the sea- 
ward rim of the shelf valleys. 

Most of the relict Pleistocene river channel network has been filled in with sedi- 
ments. However, subtle changes in relief in some areas of the shelf surface of the 
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Virginian Sea are suggestive of former channels. Examples of these transverse shelf 
valleys are found between the mouth of Chesapeake Bay and Norfolk Canyon (Susquehanna 
Valley), from the Delaware Bay shelf valley to the shelf edge (Delaware Valley), from the 
Chesapeake Bay shelf valley southeastward to the shelf edge (Virginia Beach Valley), 
from the Oregon Inlet, N.C., vicinity southeastward to the shelf edge (Albemarle Valley), 
and from the Metomkin-Assawoman Island vicinity east- southeastward to Washington 
Canyon. The valley names are adopted from reference 99. The dimensions and gradients 
of these submarine canyons (from ref. 144) are compared in table V with subaerial 
canyons. 


TABLE V.- COMPARISON OF DIMENSIONS AND GRADIENTS OF CANYONS ABOVE AND BELOW SEA LEVEL 

jpata from reference 144, pp. 122-123J 
(a) Dimensions 


Canyon 

Mapped 
total 
length, 
n. mi. 

(a) 

— 

Length 
of part 
invading 
shelf, 
n. mi. 

Maximum 
width, 
n. mi. 

Maximum 

depth 

below 

adjacent 

ridges 

Greatest 
mapped 
depth of 
floor below 
sea level 
(a) 

Mean 

Longitudinal 
gradient in 
shelf area 

Mean 

longitudinal 
gradient in 
belt of 
continental 
slope 

meters 

feet 

meters 

feet 

Hudson 

43 

15 

7 

1219 

4000 

2256 

7 400 

1:33 

1:40 

Wilmington 

56 

9 

6 

1097 

3600 

2865 

9 400 

1:18 

1:53 

Baltimore 

74 

9 

5 

823 

2700 

3048 

10 000 

1:18 

1:66 

Washington 

56 

6 

3 

732 

2400 

2499 

8 200 

1:20 

1:53 

Norfolk ....... 

52 

10 

3 

914 

3000 

2743 

9 000 

1:19 

1:48 

Colorado River . . . 

--- 

— 

13 

1859 

6100 

— 


1 : 530 (mean gradient) 

Snake River 

— 

— 

9 

2408 

7900 

— 


1 : 400 (mean gradient) 


a Actual dimension is greater. 


(b) Average longitudinal gradients of canyons (by groups) 


Regional group 

Average number 
of canyons 

Gradient of canyon as a whole 

Gradient at head of canyon 

percent 

meter s/n. mi. 

feet/mile 

percent 

meter s/n. mi. 

feet/mile 

Eastern United States . . . 

14 

5.5 

102 

290 

7.3 

135 

385 

Western United States . . . 

29 

4.83 

89 

255 

9.96 

184 

525 

Eastern Asia 

28 

7.0 

130 

370 

14.4 

267 

760 

Indian Ocean 

5 

9.8 

181 

517 

20.0 

377 

1075 

Mediterranean Sea 

12 

10.9 

202 

575 

15.2 

281 

800 

Oceanic Islands 

4 

13.8 

266 

759 

21.0 

395 

1125 

Off large rivers 

9 

1.7 

32 

90 

3.24 

60 

170 

Colorado River 

--- 

0.2 

4 

10 

— 

--- 

— 

Snake River 

— - 

0.25 

5 

13 

— 

— 
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Virginia Beach Massif .- Virginia Beach Massif, between the Susquehanna Valley and 
the Virginia Beach Valley, is an extensive shallow, relatively level- topped topographic high 
lying approximately between the depth contours of 18.3 and 21.9 meters (60 and 72 feet). 
(See fig. 1.) This imposing large-scale relict feature, of probable interfluve origin, con- 
tains a superimposed irregular ridge and swale bathymetry, which is delineated by the 
depth contour of 18.3 meters (60 feet). The Virginia Beach Valley, flanked to the north- 
east by the Virginia Beach ridges on the topographic high and to the southeast by the False 
Cape ridges, is indeed suggestive of a series of relict ebb-tidal deltas formed as the sea 
level rose and the estuary mouth retreated, as hypothesized by Swift, Kofoed, Saulsbury, 
and Sears (ref. 99). 

This complex topographic high, originating as an interfluve feature, with subsequent 
superimposed tidal-delta- associated ridges, that have been modified under the present 
shelf hydraulic regime, has been named the Virginia Beach shoal retreat massif by Swift, 
Kofoed, Saulsbury, and Sears (ref., 99); this name has been adopted for this study. 

Small-Scale Shelf Relief Elements 

Linear ridges . - Superimposed on the larger relief elements is an undulating ridge 
and swale bathymetry composed of shoals with less than 9.1 meters (30 feet) of relief, 
with the long axis generally extending from 1 to 10 miles and oriented such that they form 
a small angle (peak at 35°) with the present shoreline (ref. 145). These shoals are 
thought to have formed under the present shelf hydraulic regime because marked seismic 
and grain-size discontinuities exist between the shoals and the underlying strata which 
are generally older than 7000 years (refs. 145 and 146). Moreover, the mineralogy and 
granulometric characteristics of many of the shoals are often directly related to the 
beaches along the adjacent shoreline (ref. 145). 

Linear ridges, separated by valleys called swales (ref. 80), are most prominent 
opposite the shorelines of Delaware and Maryland, the southern Delmarva Peninsula, the 
Virginia-North Carolina state line, and Oregon Inlet to Rodanthe, N.C. 

The depth and orientation of over 200 of the linear ridges on the U.S. East Coast 
Continental Shelf is shown in figure 11 (data from ref. 145). Note the bi modal depth dis- 
tribution with clusters of shoals at depths of 6.1 to 9.1 meters (20 to 30 feet) and 12.2 to 
16.8 meters (40 to 55 feet) (and possibly a third mode at depths greater than 24.4 meters 
(80 feet)). These depths do not appear to be related to depths of prominent terraces; 
instead, they may be related to depths at which the most frequent waves begin to apprecia- 
bly interact with the sea floor. (Compare fig. 11 with figs. 5, 7, 8, and 10.) The right 
histogram in figure 11 shows the azimuth distribution of the same 200 linear ridges, with 
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Class boundaries, deg 



major axis of the shoals having a mean azimuth (i.e., compass direction) of 32°. Two 
modes are suggested at approximately 5° and 35°, with a third mode possibly at -30° 

(i.e., 330°). 

Arcuate shoals .- The arcuate shoals are most prominent when associated with 
capes such as within Chincoteague Shoals opposite the south end of Assateague Island, 

Md. They are even more extensive immediately south of the study area, within Diamond 
Shoals opposite Cape Hatteras, N.C, Arcuate shoals are also located opposite the mouths 
of nearly all the inlets along the coast of the Virginian Sea. The formation of the inlet 
shoals (i.e., ebb-tidal deltas) is related to the tidal- current-wave interaction, and they 
often have an important effect on the nearshore wave refraction patterns (ref. 25). 

Probably the largest arcuate shoal in the study area is one associated with the 
entrance to Chesapeake Bay. Though highly bisected and cut by tidal channels, the dis- 
tinct convex -seaward arcuate shape of this intermittent sand body, encompassing the 
mouth of the Bay, can be delineated from the detailed bathymetry. This huge sand body, 
suggestive of an ebb-tidal delta, may also be directly related to the origin of linear ridges 
adjacent to False Cape. Indeed, many of the linear ridges, especially those attached to 
shore, as well as many of the arcuate shoals may owe their origin, in part, to the forma- 
tion of now relict ebb-tidal deltas. 


STUDY RESULTS 


Variations in Wave Behavior With Different Combinations of Input Parameters 

Because of the very large number of combinations of wave direction, period, tide 
and height input conditions, large variations in wave behavior may result in any particular 
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portion of the shelf and along the shoreline. Also the large areal extent of the study 
adds to the complexity of the voluminous data. (A detailed analysis of these data is in 
progress.) 

However, in spite of the complexity of the computer -generated data, certain definite 
features and patterns of wave behavior are clearly displayed by the wave refraction dia- 
grams (figs. B1 to B70). One of the most notable of these patterns is a definite spatial 
periodicity in wave energy distribution over the shelf and along the shoreline. This 
periodicity results from the interaction of the ocean waves with the numerous relief ele- 
ments on the continental shelf. These shelf relief elements include the shelf-edge can- 
yons (for 12-second or longer waves) such as Washington and Norfolk Canyons; shelf 
valleys such as at the entrances to Delaware and Chesapeake Bays; ridge and swale 
bathymetry, most notably adjacent to the Delmarva Peninsula and Virginia Beach, Va.; 
and shore- connected northeast-oriented ridge systems such as Bethany Beach, Del.; 

False Cape, Va.; and Rodanthe, north of Cape Hatteras, N.C. 

The resultant shoreline wave energy distribution varies with wave approach direc- 
tion, wave period, stage of the tide, and changes in sea level from the inverse barometric 
effect associated with moderate to severe storms. This spatial variation in wave energy 
distribution (i.e., alternate zones of wave energy concentration and diminution) also varies 
directly with the period of the incoming waves. For waves of 6 to 8 seconds, these peri- 
odic zones are 1 to 5 n. mi. in length along the shoreline, and for waves of 12 to 14 sec- 
onds, these zones may be 15 to 25 n. mi. in length. The variations in the widths of these 
zones appear to be related to the distance from shore that the waves begin to interact with 
the shelf relief elements. This spatial wave energy distribution along the shoreline will 
affect the morphology and the long-term erosional history of the shoreline. Obviously, 
documentation of these trends needs to be undertaken in any shoreline planning. 

Ocean- W ave — Continental - Shelf Interaction 

Shelf- edge canyons .- Norfolk and Washington Canyons head at depths of approxi- 
mately 87.8 and 73.2 meters (288 and 240 feet), respectively. They, therefore, appear to 
have an important effect on waves with periods of 12 seconds or greater. These effects 
are twofold: 

(1) Relatively decreased concentration of wave energy downwave from the canyons 

and the creation of a shadow zone of low wave energy along the shore about 

20 n. mi. in length due to a divergence of the wave rays 

(2) Greater wave energy on either side of this shadow zone as a result of the con- 

vergence of the wave rays refracted away from the shadow zone 

These effects have been noted by Munk and Traylor (ref. 147) along the California shore- 
line where, because of the closeness of the canyons to shore, the effects are more 
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discernible. Though apparent for 12- and 14-second waves from all directions (see 
diagrams in appendix B), the effects of the canyons on the waves will vary with wave 
approach direction because of additional wave -shelf interaction in the downwave direc- 
tion. An example of this is the shoreline effect of Norfolk Canyon, which for northeast 
waves is slight because of additional refraction but much more noticeable for east- 
northeast and east waves. The shadow zone moves north along the shoreline as the wave 
approach direction moves south; the Norfolk Canyon shadow zone moves from the False 
Cape vicinity to Virginia Beach and then to the Chesapeake Bay entrance and Smith Island, 
with the wave approach direction changing from northeast to east-northeast to east, 
respectively. 

Shelf valleys .- Two of the larger shelf valleys are located adjacent to the mouths 
of Delaware and Chesapeake Bays and display a southeast orientation. The Delaware 
Shelf Valley has appreciably more relief and steeper sides than the Susquehanna Shelf 
Valley (fig. 1).' Like the shelf-edge canyons at the seaward side of the shelf, these relict 
features were probably formed by subaqueous stream processes when sea level was 
lower. Present interaction with the shelf valleys results in one of the most important 
effects of wave-shelf interaction, that is, the development of confused seas over the sea- 
ward margin of these valleys. This effect is due to landward moving waves encountering 
an abrupt increase in depth and is most dramatically illustrated for the Delaware Valley 
by 12- second waves from the south-southeast. As can be readily seen in the diagrams 
(appendix B), subtle changes in wave approach direction and wavelength will drastically 
alter the "strength” of the confused sea, which is referred to technically as a straight 
caustic. Straight caustics have been demonstrated in a wave tank and described theoreti- 
cally by Chao and Pierson (ref. 41). The concentration of wave energy from southeast 
swell (as exemplified by increased wave heights) along the seaward side of the Delaware 
Valley may have been observed and photographed by LaRC personnel in August 1973. 

Ridge and swale bathymetry .- Originally described by Uchupi (ref. 80), the detailed 
morphology and suggested processes of origin and maintenance of these ridges has been 
recently summarized by Duane, Field, Meisburger, Swift, and Williams (ref. 145). A 
mechanism for maintenance of these ridges by present wave processes has been sug- 
gested by Goldsmith (ref. 3 1), in which wave refraction around a topographic high would 
cause the sediment to be redeposited downwave from the high (fig. 12). Thus, such ridges 
would have their long axis oriented parallel to the wave rays from the dominant wave 
approach direction, as is the case for most of the ridges on the shelf. 

Irrespective of the genetic mechanism of this ridge and swale bathymetry, an 
examination of the wave refraction diagrams clearly shows that this bathymetry has the 
most significant overall effect of any of the shelf relief elements on the ocean waves. As 
a result of wave refraction over this regular ridge and swale bathymetry, alternating 
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zones of relative wave energy concentra- 
tion and diminution are produced downwave 
from this bathymetry, both over portions 
of the shelf and along the shoreline. 

Shoals, banks, and massifs .- Exten- 
sive shoal areas, forming topographic 
highs on the shallow continental shelf, will 
also greatly affect incoming wave patterns. 
Such areas are found offshore from Ocean 
City, Md. (Gull Banks and Fenwick Shoals), 
the south end of Assateague Island, Md. 
(Chincoteague Shoals), Parramore Island, 
Va. (Parramore Banks), Virginia Beach to 
False Cape, Va. (Virginia Beach Massif), 
and Oregon Inlet and Rodanthe, N.C. (Platt 
and Wimble Shoals, respectively). 



An example of the effects of these Figure 12.- Schematic illustrating a mechanism 

offshore Shoal areas on nearshore cir- proposed to explain growth and maintenance 

onsnore snoai aieas on nearsnore cir of linear ridge on shelf and its effect on 

culation patterns can be seen in the wave energy distribution along shoreline. 

vicinity of Virginia Beach, Va., which is < ‘ lrom ref. 31.) 

greatly affected by the adjacent, extensive Virginia Beach Massif. Here, the waves with 
periods of 10 seconds or shorter from the north -north east, northeast, and east-northeast 
are, for the most part, refracted away from the resort area by the Virginia Beach Massif 
to the Chesapeake Bay entrance and the Back Bay- False Cape area. In a similar manner 
waves from the east- southeast, southeast, and south- southeast are concentrated in the 
Virginia Beach and adjacent offshore area. These phenomena result in the dominant 
northward longshore drift observed in the Virginia Beach area; this might be because 
greater wave energy reaches the area from the southern quadrants than from the north, 
resulting in a net nearshore sediment transport to the north. Harrison, Brehmer, and 
Stone (ref. 148) suggested that the observed northward sediment transport in the Virginia 
Beach area was due to a large nontidal eddy related to the circulation originating at the 
mouth of the Chesapeake Bay. It should therefore be noted that both effects may be 
occurring and that neither the wave or current-induced circulation patterns are mutually 
exclusive. 


Shoreline Wave Energy Distribution 

Shore-connected linear ridges .- Of critical importance to the final shoreline wave 
energy distribution and intimately related to shoreline changes, are the several series of 
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shore- connected linear sand ridges, which, like the linear ridges offshore, generally have 
their long axes alined to the northeast. Extensive shore- connected ridge systems are 
found off Bethany Beach, Del.; False Cape, Va.; and Oregon Inlet to Cape Hatter as, N.C. 

Moody’s hand-drawn refraction diagrams for the Delaware ridges (ref. 149) indi- 
cated shoreline wave energy concentrations over the ridge crests for all common wave 
approach directions. These concentrations were apparently substantiated by Moody with 
measurements of the irregular beach changes at Bethany Beach following the March 1962 
storm. 

Vincent (ref. 150), in a very detailed statistical analysis of shoreline meanders 
along the Rodanthe-Cape Hatteras shoreline, found a definite relationship between the 
Rodanthe shoreline rhythms and the adjacent ridge and swale bathymetry, which could be 
related through the wave refraction patterns. 

Though the detailed wave refraction patterns over these ridge systems are smaller 
than the scale of these diagrams, it is significant with respect to hypotheses advanced 
concerning the origin of these shoals in reference 100 to note that the shorelines of all 
three of these ridge systems are areas of above average wave energy concentrations 
from many different directions. These concentrations are most apparent for waves with 
periods of 8 seconds from the northeast and southeast and to a lesser extent for waves 
with periods of 6 and 10 seconds from these and other directions. 

A second-order grid of the area offshore from southeast Virginia is now being 
examined in some detail at VIMS to determine the wave refraction patterns over the 
False Cape ridge system. 

Inlets and wave energy concentrations .- There are definite shoreline areas with 
pronounced wave energy concentrations for several wave approach directions. Two of 
the most obvious examples are the False Cape area and the Chesapeake Bay entrance. 
This concentration appears to be a direct result of the extensive Virginia Beach Massif 
which is relatively shallow, 18.3 meters (60 feet) deep, and close to shore (10 n. mi.). 
Present studies are following these waves into Chesapeake Bay in order to study the 
resulting refraction patterns and delineate the contribution of these ocean waves to the 
wave climate of the Bay. 

One of the more intriguing results of the present study is the close proximity of 
the coastal inlets and the shoreline areas of wave energy concentration; this suggests a 
causal relationship. This relationship is especially true for the numerous inlets of the 
southern Delmarva Peninsula (e.g. , Wachapreague Inlet) and for Oregon Inlet; thus, the 
inlets may owe their existence, in part, to wave action because the inlets are in areas of 
wave energy concentration. Also related to this is the apparent spacing of wave energy 
concentrations for waves with periods of 8 and 10 seconds for various directions, which 
is remarkably close to the spacing of the inlets. This shoreline spacing, as indicated by 
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the convergences of wave rays, is directly related both to the depth and relief of the 
various shelf elements and the length of the waves that would be most directly influenced 
by these relief elements. The relatively shallow, wide shelf helps to emphasize these 
relief elements. The slight wave refraction of waves with periods of 4 and 6 seconds 
and the abrupt increase in refraction for waves with periods of 8 seconds and longer, for 
all wave approach directions, is clearly related to this morphogeometry of the shelf. 

An excellent example of the regular periodicity that can result from the ridge and 
swale bathymetry is shown in the wave refraction diagrams (figs. B38, B20, and B26) for 
waves with periods of 10 seconds from the southeast for the North Carolina shoreline and 
from the east-northeast and east for the Eastern Shore of Virginia. Here, shoreline 
zones of wave ray divergences and convergences alternate at approximately intervals of 
5 n. mi. A detailed study has been made of such alternating zones of wave energy along 
20 n. mi. of shoreline off the southern Delmarva Peninsula (ref. 25). A noteworthy cor- 
relation was observed between observed rates of shoreline erosion since 1852 and areas 
of wave energy concentrations as computed in this second-order Virginian Sea Wave 
Climate Model (fig. 13). 

Other areas of relatively high shoreline wave energy concentrations for specific 
wave conditions are between Indian River Inlet, Del., and Ocean City, Md.; the south end 
of Assateague Island near the Virginia- Maryland state line; Cobb Island to Smith Island, 
Va.; and Corolla to Kitty Hawk, N.C. 

Thus, shoreline wave energy concentrations are formed by wave refraction over 
the various shelf relief elements. The location and spacing of these shoreline zones of 
wave ray convergences and divergences will depend on 

(1) The wave approach angle 

(2) The size and extent of the adjacent shelf relief elements 

(3) The distance from shore and depth of these relief elements 

(4) The regularity and rhythmic nature of the relief elements 

(5) The interaction of wave refraction patterns from more than one set of relief 

elements 


Prediction of Shoreline Changes 

Figure 13 illustrates the surprisingly strong correlation near Wachapreague Inlet, 
Va. (between Cedar and Parramore Islands) between shoreline variations in measured 
rates of erosion between 1852 and 1934 and areas of wave energy concentrations in 1852 
determined by using 1852 bathymetry as input data. The three wave energy peaks for the 
1852 erosional waves with periods of 4 seconds show a strong correlation in both location 
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Figure 13.- Historical shoreline changes along the Eastern Shore of Virginia compared 
with changes in shoreline wave energy distributions. 1852 to 1934. 


and magnitude with the three areas of shoreline erosion. Also, the three wave energy 
peaks for the 1934 accretional waves with periods of 12 seconds show a similar correla- 
tion with the three areas of shoreline accretion. In each of these two cases, the amounts 
of erosion are also directly related to the magnitude of the computed wave energy. Coin- 
cident with the irregular shoreline changes that resulted in the increased development of 
the offset inlets, the shoreline wave energy for the waves with periods of 4 seconds tended 
to become more evenly distributed along the shoreline between 1852 and 1934. Therefore, 
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one might be emboldened to predict that when the alongshore shoreline wave energy dis- 
tribution becomes evenly distributed, the inlet offsets will cease growing and rates of 
shoreline erosion and accretion will become more evenly distributed along the shore. 

Thus, by this evaluation of historical wave conditions based upon historical shore- 
line and bathymetric data, changing shoreline wave energy distributions with time (as a 
reflection of offshore bathymetric changes) can be delineated, and predictions can then 
be made of future shoreline and nearshore changes. Therefore, the usefulness of the 
VIMS-LaRC Wave Climate Model is not limited to present conditions, but the model is a 
valuable tool which can also be utilized in understanding the past in order to predict the 
future wave energy distribution over the shelf and along the shoreline. 

Tidal Effects 

In order to test the effect of changing water level on these refraction patterns, 

1.2 meters (4 feet) 4 were added to each of the depths and the waves from the northeast, 
east, and southeast were recomputed. 

In general, for the east and southeast directions the wave rays are less concen- 
trated along the shore for high -tide conditions, though the areas of wave energy concen- 
tration are not significantly different for equivalent low- and high-tide conditions. For 
waves from the northeast the situation is more complex, with some significant changes 
between low- and high-tide conditions involving focusing of wave ray convergences and 
changes in the locations of this focusing of up to 2 or 3 n. mi. along the shore. 

APPLICATIONS 

An important aspect of this applied research supported in part by the NOAA Sea 
Grant Program is the concomitant development of applications for the new research 
technology. 

The presentation arid analysis of the data derived from the VIMS-LaRC Wave 
Climate Model further substantiate the general recognition that ocean wave refraction 
plays a dominant role in controlling the distribution of wave energy over the continental 
shelf and along our coasts; furthermore, these data illustrate how knowledge of a specific 
refracted wave input at a site is critical to the successful implementation of shoreline 
defense programs, knowledge of continental shelf sedimentation processes, and construc- 
tion and maintenance of coastal and offshore structures. Thus, a number of agencies 
have interests in the continental shelf waters (e.g,, VIMS, NASA, U.S. Army Corps of 
Engineers, NOAA, Environmental Protection Agency, Geological Survey, Fish and Wildlife 
Service, Virginia Marine Resources Commission, Commonwealth Division of Planning and 

^Equivalent to a spring tide over most of the area. 
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Community Affairs, and Virginia Highway Department) and are potential users of informa- 
tion from analytical wave models such as the VIMS-LaRC Wave Climate Model. Also, 
industries which may benefit from wave predictions on the shelf include shipping, oil, 
fishing, and recreation. 

Although the present model is limited to a first-order linear approach, it is an 
invaluable tool for the scientific understanding of a number of shelf processes. Addi- 
tional experimental verification of this model (or an improved version) offers the oppor- 
tunity for a variety of additional potential uses. 

Advances in analytical modeling may be expected to interact with technology 
advances from other areas of science. Although all potential uses cannot be defined, 
major categories of specific long-term applications discussed in this section are 
(1) advanced coastal wave forecasting, (2) improved environmental quality, and (3) data 
supply for government, industry, and other scientific disciplines. 

Advanced Coastal Wave Forecasting 

Considerable potential exists for the application of continental shelf wave refrac- 
tion modeling for the improvement of marine forecasts of wave heights over the continen- 
tal shelf. Two major approaches to this problem are discussed: (1) the continental shelf 
modeling may be incorporated into existing Weather Service forecast techniques or (2) it 
could become part of the data analysis procedure for future satellite systems. 

The present NOAA marine forecasting technique (ref. 151) utilizes the National 
Meteorological Center (NMC) 1977-point grid system for the prediction of wind -gene rated 
waves and swell. Computations are made at each grid point spaced 180 n. mi. apart. 

Wave forecasting techniques applicable to the deep ocean are used in these computations. 
Wave-height contour charts are then derived based on values at each grid point over the 
entire North Atlantic Ocean. Because of the large grid spacing and utilization of deep- 
ocean computational methods, the present forecast system has limited application to the 
continental shelf regions where bottom topography modifies ocean swell characteristics. 

It is reasonable to expect that wave refraction models of small areas of the continental 
shelf could be incorporated within the large-scale NMC grid squares in coastal regions. 
Output (height, period, and direction) from the NMC grid could be smoothed and input to 
the refraction model which in turn would predict height distribution in finer detail over 
the continental shelf region. This same procedure is essentially followed for a smaller 
geographic region in reference 42. 

An alternative approach to more accurate wave forecasts would be to utilize a sys- 
tem similar to the prospective Seasat oceanography satellite (ref. 152) coupled to the 
wave refraction model by using microwave instrumentation for all-weather capability. 
Deep-ocean wave directional spectrum data could be input to the refraction model for 
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prediction of continental shelf swell conditions. Wave spectra and/or wind data could be 
used to predict wind- generated wave height distribution. Recent successes in monitoring 
wave and wind characteristics with microwave instrumentation (refs. 153 to 155) give 
encouragement that such advances are possible. Harrison and Green (ref. 156) show that 
satellite orbits are possible which maximize oceanographic and wind observations in the 
coastal zone. 

As the continental shelf regions of the United States are subject to increased utili- 
zation, there will be a need for more detailed wave forecasts over the region just as 
detailed meteorological forecasts are now required around airports. A study by Pierson 
in 1972 (ref. 23) has shown that wave refraction effects over high-relief bathymetry in 
the North Sea could have been directly related to the sinking of two fishing boats. A sim- 
ilar experience may have been encountered by Florida fishermen off Cape Charles, Va., 
in December 1973 (ref. 157). Coping with the natural weather and sea environment is 
probably the largest single problem in offshore operations (ref. 158). It is clear that as 
advanced marine prediction systems are developed, large-scale analytical wave refrac- 
tion modeling will play an increasingly important role. 

Improved Environmental Quality 

Although prediction of surface wave conditions is of obvious value to offshore oper- 
ations, this same capability has important impact on the ability to analyze a number of 
environmental problems. The ability to analyze and predict shoreline changes caused by 
either natural or manmade bathymetry changes on the shelf has been discussed in the 
section "Study Results." Prediction of surface wave conditions also has potential impact 
on ocean dumping operations. For products categorized as liquid wastes, it has been 
suggested that disposal sites be selected which allow rapid dispersal near the surface. 
This implies that it may be desirable to dump these products in a region with high sea 
state where turbulent mixing is high. Thus prediction of weather and wave conditions 
could be used to vector dumping vessels to high energy regions of the shelf. 

A theoretical study utilizing second-order Stokes wave theory (ref. 159) suggests 
that wave-induced currents may, under some conditions, be the same order of magnitude 
as wind-driven currents on the Mid-Atlantic Continental Shelf. This suggestion has 
important implications concerning the transport of pollutants both on the surface and at 
depth in shelf waters. Predictive capability of the bottom currents is particularly impor- 
tant to site selection for disposal of both inert and degradable solid wastes. Surface cur- 
rent conditions are important for prediction of oil-spill movement. Close to land, wave- 
induced longshore currents are important to the transport of effluents of sewage outfalls. 
One method for monitoring shelf circulation (noted in ref. 159) may be the application of 
wind, wave, and Gulf Stream data from remote sensors to analytical circulation models 



which would calculate both surface and subsurface currents. Such modeling should 
include wave refraction calculations as well as geostrophic balance, wind stress, salt 
balance, thermal balance, and fresh -water runoff. 


Data Supply for Government, Industry, and 
Other Scientific Disciplines 

As frequent, all-weather wave observations become available from systems such as 
the planned Seasat satellite, a greater scientific understanding of wave conditions on the 
continental shelf will evolve. This understanding will enable more sophisticated analyt- 
ical models. For example, it may be possible one day to transform wave directional 
spectra statistics taken at one location to another location through a series of verified 
refraction calculations. Such operations could reduce the need for wave gage installation 
and long duration measurement periods at sites of proposed offshore construction (such 
as is now underway on the New Jersey coast). 

With respect to direct application of the Virginian Sea Wave Climate Model to oil- 
spill problems, recent studies indicate that the spread of oil on the sea surface is a com- 
plex problem, with the oil movement mainly a function of the direction and velocity of the 
surface winds and tidal currents (ref. 160). Sonu, Murray, and Smith suggested in refer- 
ence 161 on pages 17 and 18 that wave action was an important factor in the spread of oil 
under sea conditions inasmuch as wave action affected mixing in the water layer which in 
turn influenced vertical eddy viscosity and, hence, the vertical velocity profile. There- 
fore, although this model may be somewhat limited in its ability to predict the movement 
of oil, it could be combined with other models (discussed previously) and contribute basic 
data on swell conditions and areas of potential confused seas, especially in crisis 
situations. 

If data derived from either the present VIMS-LaRC Wave Climate Model or more 
sophisticated versions are used, more rapid definition of engineering specifications, 
appropriate legislation, insurance acturial data, and environmental impact effects may be 
possible. Such information is necessary for the correct management of environmental 
resources. For example, accurate knowledge of shoreline wave energy distributions is a 
must in defining coastal setback lines such as that legislated in 1971 in Florida (ref. 162). 

By using deep-water wave conditions from Marsden square or hindcasted storms, 
data from the VIMS-LaRC Wave Climate Model for specific wave conditions can be put 
into a statistical representation of shoreline wave energy distribution. Risk factors can 
then be developed by combining these data with probabilistic descriptions of storms (i.e., 
frequency and intensity). Such data should be of great interest to insurance companies 
and aid in the formulation of setback lines. 
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These wave data may be used specifically to assist in the analysis of the distribu- 
tion of the sediments of the sea floor of the Virginian Sea. In a detailed study conducted 
by Nichols (ref. 104), these sediments displayed large grain- size variations on the shelf. 
Comparison of a detailed grain- size distribution map prepared by Nichols, with a shelf 
wave energy and bottom wave scour ’'map,'’ contoured for a variety of wave conditions, 
should provide valuable information as to the causes of the grain-size distribution (i.e., 
relict processes and/or modern hydraulic processes). If a comparison suggests a rela- 
tionship between the grain-size and wave energy distributions, then the wave data can be 
used to suggest areas for future sediment sampling in order to further test the hypothesis. 
Thus, it is clear that regional- size banks of data concerning wave orbital velocities and 
wave-induced currents will provide beneficial information to scientists in other 
disciplines. 


CONCLUDING REMARKS 

A Wave Climate Model for the Mid- Atlantic Continental Shelf and shoreline encom- 
passing 20 000 square n. mi. and 200 n. mi. of shoreline has been computed for 122 dis- 
tinct wave and tidal conditions. Computations of 19 wave parameters along more than 
30 000 wave rays have resulted in the generation of over one million pieces of useful data, 
which will be used as 

(a) A data bank of useful wave information available at the Virginia Institute of 

Marine Science (Gloucester Point, Va. 23062) to user agencies, other officials, 
and companies involved with the development and use of the continental shelf 
and shoreline 

(b) A source for additional, more detailed analyses, such as testing of proposed 

modifications prior to initiation of projects 

(c) Input to smaller closer spaced near-shore grids 

(d) The nucleus for an early warning system from severe storm events for coastal 

areas 

The most important wave process on the continental shelf is the interaction between 
the ocean waves and the various shelf relief elements. This process results in nonuni- 
form wave energy distributions over the shelf and along the shoreline. Detailed analyses 
of these ocean- wave— -shelf interactions for specific shelf and shoreline areas should 
provide valuable input for the understanding of the present processes, the past develop- 
ment, and future changes of the continental shelf and coastal geomorphology. 
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SUMMARY OF MAP PROJECTION CONSTRUCTION 

By Norman T. Carpenter 
Clark University Cartography Laboratory 
Worcester, Massachusetts 

Construction of the Map Projection 

For the construction of the map projection, a grid of 0.0 508- meter (2-inch) squares 
was scribed on Keuffel & Esser 0.127-millimeter (0.005-inch) scribing film, 2.74 meters 
(108 inches) long and 1.016 meters (40 inches) wide. The measurements were made by 
surveyor’s metal tape for intervals of 0.3048 meter (1 foot) or more and by draftsman's 
scale for shorter intervals. 

The center intersection on the grid was designated the origin of a rectangular coor- 
dinate system, with the ordinate parallel to the long dimension of the grid. The origin 
was designated lat. 37° N, long. 75° W in the geographic grid. A computer program 
(ref. 62) calculated abscissa and ordinate values in inches for intersections at 10- minute 
intervals of latitude and longitude, based on the sphere of equal area, radius equals 
6 370 997 meters (3440 n. mi.). The intersections were plotted by hand and connected 
by scribed interrupted lines. The estimated accuracy of location of all intersections is 
within 0.508 millimeter (0.02 inch). 

A positive contact print made from the original scribe coat was measured. The 
grid on this print, made on a stable-base material, is 2.75 meters (108.1 inches) long 
and 1.0168 meters (40.03 inches) wide, an expansion of about 1/10 of 1 percent. The 
diagonals of the grid were equal: no deviation from rectangularity could be measured. 

The nominal scale of the projection is 2.54 centimeters (1 inch) equals 2.5 n. mi . , ^ 
or 1 : 180000. Variations in this scale are discussed in the following section. 

The map extends from latitude 35° to 39° N and from longitude 74° to 76° W. With 
no change in the present map, the projection can be extended in any direction. 

Scale Factor 

Scale factor is the term used to define the variation in areal exaggeration of a pro- 
jection. Since the scale of the projection is 1 : 180000, it corresponds to a globe 1/180000 
the size of the sphere of equal area to the Earth. All distances along the standard line 

5 For the map projection construction, 1 n. mi. = 1852 meters (6076 feet); for the 
example in this discussion, 1 n. mi. = 1829 meters (6000 feet). 
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are exactly the same on both globe and map. Leaving the standard line, distances are 
slightly greater on the map than on the globe, so that in the actual scale, the denominator 
is slightly less than 180 000. The scale factor is equal to the denominator of the nominal 
scale, or 180 000, divided by the denominator of the actual scale. The difference in scale 
is not great on this map. On the standard line, 75°, the scale factor is 1.0000 and 
2.54 centimeters (1 inch) equals 4572 meters (15 000 feet); 50 n. mi. east or west from 
the standard line, the scale factor is 1.0001 and 2.54 centimeters (1 inch) equals 
4571.5 meters (14 998.5 feet). (See fig. 1 of ref. 55.) 

Convergence 

A projection is conformal when the angular relationships are the same at any point 
on the map as at the corresponding point on the Earth; therefore, conformal maps are 
used in navigation and surveying. The price of conformality is area exaggeration, but 
that is no problem in the limited extent of the UTM. However, it should be understood 
that while all directions are true at a point on the UTM, they vary from one point to the 
next. There are two characteristics of the projection to consider in connection with 
direction and the measurement of azimuth. 

First, on this map the latitude and longitude lines have been superimposed on a 
square grid, which is more convenient to use, since locations can be expressed in 
Cartesian coordinates and directions expressed in plane trigonometry. The grid direc- 
tion parallel to the central meridian is called grid north GN. However, grid north is the 
same as true north only on the central meridian and on the equator. The angle between 
grid north and true north is called the convergence, and it increases with distance from 
the central meridian and with distance from the equator. On this map it amounts to less 
than 1°. 

Second, except for the central meridian and great circles at right angles to it, a 
straight line on the Earth’s surface is not a straight line when projected on the UTM, and 
a straight line on the UTM does not pass through the same intermediate points as a great 
circle between the corresponding end points on Earth. This difference is greatest for 
north- south lines, increasing with distance from the central meridian, but it is slight on 
this projection. Take the two points 35° N, 74° W and 39° N, 74° W. A straight line on 
the Earth’s surface between these points is, of course, the meridian 74° itself. The 
meridian is a curve on the map and at its midpoint,- 37° N, it deviates from a straight 
line by 0.302 millimeter (0.0119 inch) or 54.5 meters (178.5 feet) on the Earth. Also, 
a straight line drawn on the map tangent to the meridian at 39° N, 74° W would cross 
35° N about 1.21 millimeters (0.0478 inch) from the meridian (218.5 meters (717 feet) 
on the Earth’s surface). 
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WAVE REFRACTION DIAGRAMS 


Included in this appendix are 70 wave refraction diagrams (figs. B1 to B70) com- 
puted for the 122 wave input conditions listed in table I. Because linear wave theory was 
used, computations for initial wave heights of 0.6 and 1.8 meters (2 and 6 feet), with 
otherwise identical input conditions, will give identical refraction diagrams; therefore, 
the initial wave height was 0.6 meter (2 feet) for all diagrams shown. (To convert val- 
ues of depth contour from feet to meters, multiply by 0.3048.) The conditions for the 
70 refraction diagrams are given in the following index: 
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Figure B3.- Wave rays computed with following input conditions: 
AZ = 0°; T = 8 sec; Tide = 0. 
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Figure Bll.~ Wave rays computed with following input conditions 
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Figure B30 . - Wave rays computed with following input conditions: 
AZ = 112.5°; T = 6 sec; Tide = 0. 
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Figure B37 - - Wave rays computed with following input conditions 
AZ = 135°; T = 8 sec; Tide = 0. 
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Figure E44.- Wave rays computed with following input conditions: 
AZ = 157.5°; T = 10 sec; Tide = 0. 
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Figure B69-- Wave rays computed with following input conditions 
AZ = 135°; T = 12 sec; Tide * 1.2 m (4 ft). 
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